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We consider the single-handed spinor field in interaction with its own gravitational field described 
by the set of field equations given by Weyl field equations written in terms of derivatives that are 
covariant with respect to the gravitational connection plus Einstein field equations soured with 
the energy tensor of the spinor: for the Weyl spinor and the ensuing spacetime of Weyl-Lewis- 
Papapetrou structure, we will find all exact solutions. The obtained solution for the metric tensor 
is that of a PP-wave spacetime while the spinor field is a flag-dipole. 
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I. INTRODUCTION 

Not long after Einstein wrote the field equations of gravitation, Schwarzschild applied them outside a stationary and 
spherically symmetric matter distribution obtaining the corresponding exact solution that describes the surrounding 
gravitational field; soon after Maxwell wrote the field equations of electrodynamics, Maxwell himself applied them in 
the vacuum finding exact solutions describing the propagation of electromagnetic waves; and it took even shorter for 
Dirac to write the spinorial field equations finding in the free case the exact solutions representing the material wave 
function. All these instances refer to cases of paramount importance in physics because the knowledge of those exact 
solutions provided a detailed description of a given physical system, whether it was planetary precession and light 
bending, light-rays propagation, or the properties of electrons; on the other hand, such solutions are obtained in the 
vacuum or for free fields, that is for isolated systems. Interacting systems can of course be studied, but in that case 
the quest to find the exact solutions immediately proves to be a much tougher task to accomplish. 

A first example of such a kind of interacting system can be found again when Einstein equations are sourced by 
the energy density of an electrodynamic field whose Maxwell equations are written with the covariant derivative 
containing the spacetime connection, and for this system it was Kerr who studied the axially symmetric case finding 
exact solutions; to our knowledge however, this first example is also the last. And to an attentive examination, this 
is not a fully coupled system either, because if the gravitational and electrodynamic fields are coupled, nevertheless 
the material charged distribution is thought to be circumscribed within the singularity of the Kerr metric and it does 
not appear as the current density of the electrodynamic field equations nor as the energy density of the gravitational 
field equations. Interacting systems of matter and its own interactions can be discussed, but in this case the search 
for exact solutions has clearly become harder and harder to achieve. 

If one wishes to study the complete system of Einstein equations sourced by the energy density of electrodynamic 
and material fields, with Maxwell equations having the covariant derivatives containing the spacetime connection and 
sourced by the current density of the matter field, and the matter field equations having all covariant derivatives 
containing both the spacetime and the gauge connection, then the system becomes considerably more difficult, and 
in fact no exact solution has been found yet. 

On the other hand, looking for exact solutions of the complete and fully coupled system of matter and its own 
interactions is too fundamental a problem to be abandoned because of its difficulty. So if one were to embark on the 
enterprise of looking for such exact solutions, the search may start by considering the simplest system of all. 

The simplest system one may think to study would have to consider the spinor field in interaction with its own 
gravitational field because no field can exist without an energy density; however, considering spinors without inter¬ 
actions of electrodynamic essence is a way to simplify the problem since we do know spinors that have no electric 
charge: therefore, a first assumption we will consider is that of focusing on neutral spinor fields. 

Another possible simplification may come from the fact that spinors, despite being fundamental, nevertheless are 
not irreducible objects, but they are constituted by two chiral parts which, despite being in interaction, nevertheless 
are independent; however, their interaction is encoded in the presence of the mass term and as a consequence the 
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assumption of masslessness would allow for the complete splitting of left-handed and right-handed projections, and 
the possibility to study them separately: hence, a second assumption will be that of masslessness and the consequent 
selection of a single-handed component for the spinor field. 

In recent times, extensions of the standard model of particle physics have had to deal with the inclusion of a 
possible mechanism of neutrino oscillations, and the subsequent connection to neutrino masses, but there is a problem 
in such a scenario: although it has been demonstrated that neutrino masses would imply neutrino oscillations, the 
converse has never been proven; as such, we might be open to the possibility that neutrinos oscillate even if massless 
because of the action of some other mechanism. This possibility is made even more cogent by the fact that apart 
from the phenomenon of oscillation there is no other theoretical reason for neutrinos to be massive; further still, no 
experimental evidence of neutrino masses is present. If this scenario were true, then neutrinos would be massless, and 
they would constitute the physical prototype of the exact solutions we are searching in this paper. 

If this scenario were not true, and neutrinos were massive after all, then there would be no particle that can be 
modelled by a possible exact solution at the moment; but there could have been in the earliest stages of cosmic 
evolution. In fact, before cosmic inflation, the universe was supposed to be in a configuration of unbroken electroweak 
symmetry, where all particles were massless. Electroweak symmetry would impose an SU(2) x U( 1) gauge interaction 
but still there may be situations in which the spacetime curvature would be the dominant contribution, and there the 
particle would be modelled by the exact solution we plan to find. 

And even if there were no situation in which the spacetime curvature dominate, an exact solution for a single-handed 
spinor may still be of mathematical interest because it could be taken as the least-order contribution of the small-mass 
perturbative solution for a massive spinor field. That is, the solution for the massive spinor field could be obtained 
as a series expansion in the mass, with the exact solution in the massless case as the leading order, and solving the 
Dirac equation at the first order in the small mass term. 

As it stands, there may be several applications for a solution that, being exact, would be important also in itself. 

In this paper, we will look for such exact solution. 


II. GEOMETRIC AND MATERIAL FIELD EQUATIONS 


To begin, we start with the fundamental kinematic definitions and the form of the dynamical action we will employ. 

Throughout the paper, partial derivatives of a given function /( x h ) are denoted by f x h := ~§^k- Latin and Greek 
indices run from 1 to 4. The metric tensor of the spacetime is denoted by gij while the tetrad field associated with a 
given metric is indicated by in such a way that r= gije*e{, with = diag(—1, —1, —1,1) being the Minkowski 
metric, and with inverse ej, verifying therefore e^e^ 1 = 5j and efe^ = S in terms of Kronecker delta. Dirac matrices 
are indicated by q M and T* := e^q^ with q 5 = *q 4 q 1 q 2 q 3 and chiral representation is used. The spinorial-covariant 
derivatives of a Dirac field ^ are expressed as 

Di4> = ip x i - ttiip ( 1 ) 


where the spinorial-connection f \ is given by 


O — M pi pV-rhyk 
■— t9jh UJ i v e n e k l 1 


( 2 ) 


with uj i 11 v being the coefficients of the spin-connection associated through the relation 

T. h — u P p h e u 4- 


( 3 ) 


to the linear connection Y' ik as usual. Natural units h= c = 8nG = 1 are used. 

We will consider a Dirac field coupled to gravity in the Einstein-Sciama-Kibble theory. The Dirac Lagrangian is 


Ld = - - DiipT z tp) - mipip 


( 4 ) 


where m is the mass of the spinor: the total Lagrangian can be varied with respect to all fields yielding the Einstein 
field equations coupling the Ricci curvature R t] to the energy density 


•“Hi r, R(Ji3 — Ajj 


( 5 ) 


the Sciama-Kibble field equations coupling the Cartan torsion T^ • h to the spin density 


rji h _ c r 

1 ij ~ °ij 


( 6 ) 
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where 


Hj ■= | (ipTiDjij) - Djil’Titp) - -Ld Qij 


(7) 


and 


Si] h :=±-^{ rMTi.r,-]}^ (8) 

are respectively the energy and the spin density tensors and 

iT h D h ip + l -T h r V - = 0 (9) 

are the Dirac field equations determining the dynamical behaviour of the spinor field. 

After separating all curvatures and covariant derivatives into purely metric curvatures and covariant plus torsional 
terms, using the torsion-spin coupling equations © it is possible to substitute all torsion contributions into spin-spin 
contact potentials 0 , 0 : the significant part of the Einstein field equations results to be the symmetric one, which 
can be written in the final form as a relationship linking the purely metric Ricci curvature tensor and scalar Rjj and 
R to the purely metric energy density tensor as 

Rij - T^Rgij = s ij + ^■(V-7 5 7 t V’)(V-7 5 7 rip)gij (10) 

where the purely metric energy density tensor is 

:= ^ - D(ji> FqV-) 

and the Dirac field equations handled in a similar way yield the final expression 

iT h Dhi/j -— [(V-VO + 7 5 V07 5 ] V- ~ TO V- = 0 

16 

in terms of the Di which denotes the purely metric spinorial-covariant derivative. 

This theory is general, but of course it is possible to apply it to some specific case: if the spinor is single-handed 
then it verifies relationships of the type 7 = ±iIj and the above system of field equations m■ CU) and m becomes 

Rij = ^ (v-r - D^Ti^ (13) 

and 

iT h D h 4> = 0 (14) 


( 11 ) 

( 12 ) 


showing that for single-handed particle no torsion is present after all. 


III. SPIN EIGENSTATES IN WEYL-LEWIS-PAPAPETROU SPACETIME 


In the following, the single-handed spinor will be chosen to be left-handed (obviously, right-handed spinors would 
have the same dynamical properties), and therefore the spinor will be chosen to have the general form 

V-=(q) (15) 

where L is the left-handed spinor (the right-handed would have been the lower component). 

The two components of the semi-spinor are the two opposite spin eigenstates, so that by writing the semi-spinor as 



( 16 ) 
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the solution of the eigenstate equation 5 O 3 L = ±L assign to the upper component a the meaning of the spin-up 
projection and to the lower component 6 the meaning of the spin-down projection: eigenstates can be chosen by 
vanishing one of the projections, so that it is with no loss of generality that we can assume 6 = 0 and the remaining 
degree of freedom will be described in general by a = Se+ where 2 and (3 are real functions. 

The form of the spinor is therefore 



(17) 


and from it we may compute the 5 spinorial bilinear fields: it turns out that only the vector bilinear fields 
and t/>y^y 5 ^i are different from zero, and in particular they only have the temporal and one spatial component; what 
this means is that for single-handed spin eigenstates of spinors the spatial components of the bilinear vectors select 
a privileged direction in space, identifiable with the direction of propagation, and consequently it becomes possible 
without loss of generality to restrict the spacetime to that displaying axial symmetry. 

Such a spacetime is long known to be the Weyl-Lewis-Papapetrou spacetime, in which the metric has line element 


ds 2 = —B 2 (r 2 d6 2 + dr 2 ) - A 2 (-W dt + dip) 2 + C 2 dt 2 

where all functions A(r,9), B(r,9 ), C(r,6) and W(r, 6) depend on the r and 9 variables only. 
The spinorial-connection coefficients are expressed as 




AW r y 3 y 4 5 0 yV 


= — 


4 C 

B + rB r 
2B 


2 rB 


1 2 , AWg 3 4 
7 7 + -^-77 


(18) 

(19a) 

(19b) 


A~W r i _4 A r 13 | A 2 Wg 2 4 
v ~ ABC 7 7 2 B 1 7 + ArBC 1 7 


Aa 


2 3 

7 7 


_ ( AW r WA r 

—h 


+ 


V AB ' 2 B 
C r A 2 WW, 

2B ABC 


2 rB 

T v+f+5+j 


7 1 7 4 + 


V 4 rB ' 2rB ) 7 7 

C e A 2 WWg 


2 rB 


ArBC 


2 4 

7 7 


(19c) 


with 1 l) x i ^ 0 since for spinors a possible dependence of the azimuthal angle and the time may always be 
With this form of the metric, we can evaluate the Einstein equations (fl3ll getting 


A 2 W 2 
2 C 2 


^4 j‘ /i, 


^41? 


CB 


A rr 

~A~ 


C rr 

~c 


B 


rr 

B~ 


B? 

B 2 


B r AgBg CgBg 

rB r 2 AB r 2 CB 

Bee B 2 
r 2 B r 2 B 2 


(19d) 

present. 

(20a) 


A 2 W$ 
2 C 2 


rC r 

C 


Cm 

C 


B 


rB r 

~B B BC 
rA r Agg r 2 
A 


BgCg , r 2 B 2 


+^i- 

B 2 A 


BC B 2 

AgBg 


AB 


+ 


AB 


= 0 


(20b) 



A 4 W 2 

A 4 Wg 

AA r C r 

AA rr AA r 

AAgCg 

A Agg Az? (iy 

n 


2 B 2 C 2 

2 r 2 B 2 C 2 

B 2 C 

B 2 rB 2 

r 2 B 2 C 

r 2 B 2 2 



A 4 W 2 W 2 

A 4 W 2 W% 

AW 2 A r C r 

AW 2 A rr 

3 AWA r W r 

A 2 WW r C r 

A 2 WW rr 


2 B 2 C 2 

2 r 2 B 2 C 2 

B 2 C 

B 2 

B 2 

+ B 2 C 

B 2 

A 2 W 2 

AW 2 A r 

A 2 WW r 

AW 2 AgCg 

AW 2 Agg 

3 AWWgAg 

A 2 WWgCg 

A 2 WWgg 

2 B 2 

rB 2 

rB 2 

r 2 B 2 C 

r 2 B 2 

r 2 B 2 

+ r 2 B 2 C 

r 2 B 2 


A 2 Wg ( CC r A r CC rr 


2 r 2 B 2 


AB 2 


B 2 


CCr CCgAg CCm _ AWE 2 Pt CE 2 /3 t = 0 


rB 2 r 2 AB 2 


iB 2 


(20c) 


(20d) 


2 


2 
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2 C 2 


A g B r 

■ , A r B g 

B r C e 

A r e 

Cre 

AB 

AB 

CB 

A 

' C ' 


A 2 W e E 2 

AB e E 2 

A e E 2 , Al 


16 rC 

8rB 

8r ' 


+ —r + 


CB rA rC 


= 0 


= 0 


l 2 C/3 r E 2 (4 AW(3 r r - AW g - 2 WA g +2C g ) | A 2 WW g E 2 + S 2 ( CB g - AWB e ) __ Q 


16 r 


16rC 


8 rB 


rAB r z? rA r E 2 rA 2 W r E? AS 2 
-'- 1 -'- 1 - ' - 1 - 


8 B 


16 C 


A 


~2 


= 0 


rC r E 2 + (—rCE 2 + rAE 2 W) B r rWA r z? (rACE 2 +rA 2 WE 2 )W r S 2 (— Af3 g + 2) C 


8 B 


16C 16 

S 2 {4W/3 e - 2W) A 

16 


(20e) 
(2 Of) 

(20g) 

(20h) 

(20i) 


= 0 


E 2 Cl3 ip AWE 2 /3p AE 2 (3 t A 4 WW 2 A 4 WW$ AWA r C r AWA rr 3 AA r W r 


+ 


2 B 2 C 2 + 2 r 2 B 2 C 2 


B 2 C 


B 2 


2 B 2 


A 2 C r W r , A 2 W rr AWA r : A 2 W r , AWA e C e , AWA gg , 3 AA g W g A 2 C g W g 


2B 2 C 


2 B 2 


rB 2 


2rB 2 


r 2 B 2 C 


r 2j^2 


2 r 2 B 2 2 r 2 B 2 C 
, A 2 W gg 


(20j) 


2 r 2 B 2 


= 0 


while with the spinorial-connection it is possible to expand the spinorial derivative so that Dirac equations (1141) become 

* WE 


A 


C 


C 


= 0 


\ Sft WEP V 


A 

EB r 


C 


C 

jA r 


= 0 

AW r E EC r 


rB 

E/3r 

B 


B 


2 B 2 2 rB 2AB ABC 2 BC 

9 . ^A g AEW g EC g 


= 0 


rB + 2rB 2 ' 2 rAB ArBC ' 2 rBC 


= 0 


(21a) 

(21b) 

(21c) 

(21d) 


where the separation between real and imaginary parts has already been done. 

The first four of (EU1) are the diagonal part of Einstein equations and they are those expected to determine the 
dynamics, while the remaining equations represent the off-diagonal part and they are expected to furnish constraints 
on the various components: in fact, there are 6 independent functions of the spacetime point and 14 independent 
equations, so that it is clear that some constraints on the functions must eventually be found. 

In order to solve equations m and we start by deriving (3 r and j3 g from equations (I20fl) and (I20hl) getting 


AW g A g B g 


f3 r — , „ + ~—7 — 


4 rC 

rAW., 
4 C 


2 rA 2 rB 


rA r 

~2A 


rB r 

~2B 


(22a) 

(22b) 


which can be inserted into (121cl) multiplied by BE to yield equation 


2 C 


= 0 


A , 
~A 


(23) 
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from which we get the relation 


_ = Fi(0,ip,t) 

aVc 


(24) 


in terms of F\ as an arbitrary function independent on r; also, inserting (122a|) . (I22bll and (l24l) into (I21d|> . we obtain 


(Fi) e 

rBAs/C 


= 0 


(25) 


from which it follows that F\ is only function of p and t : all in all, we obtain 

„ = F/p,t) 

As/C 


(26) 


showing that in the spinor module the dependence on azimuth and time and that on radial coordinate and eleva¬ 
tion angle are encoded within two factorized functions. Inserting all these results into equation (121 all , after some 
calculations, we get the relation 


C(r,0) 

A(r,0) 


W (r, 6) 


(Fi)t 

(Fi) v 


(27) 


from which we have necessarily 


W = A 
A 


(28a) 


(F/t _ ' 

(Fi) v 

for some constant A; equation (|28bl) implies F\ (p, t ) = F 2 (p + A t) for some arbitrary function F 2 and thus 

C = AW+ XA 


so that 

_ _ F 2 (v? + A t) 
As/C 


(28b) 


(29) 


(30) 


in which the spinor module the dependence on azimuth and time has been fixed. Now, inserting the content of 
equations (122a[> . (122bl) . (l29ll and (l30l) into (121bl) we get the further equation 


which gives 


F 2 (ip + At) — /3 1 ) 

A(A{W + A))5 


Pt = xp v 


(31) 


(32) 


fixing the interdependence of azimuth and time also in the phase of the spinor field. Collecting all results obtained so 
far, we have that if (122aD , (I22bl) . Kill) . (1301) and (15^1) hold then all Dirac equations (1^11) are satisfied automatically. 

Moreover, Einstein equations (I20fl) . ( |20g[), (120hl) and (I20il) are satisfied automatically as well. In addition to this, 
substituting (|22aj) . (I22bf> . (l29ll . (l30ll and (13211 into the remaining Einstein equations, differentiating with respect to p 
and multiplying all by 2AC/F 2 , we obtain a unique equation given by 

F 2 Pp V + 2/3^2 = 0 (33) 

where F." 2 denotes the derivative of F 2 : this equation is integrated as 


F 3 (r, 9, t ) 
F 2 


V — 


(34) 
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with F 3 (r, 9 , t) being an arbitrary function. Inserting (1341) into Einstein equations and performing the derivative of 
the so obtained equations with respect to t we can deduce that ( F 3 ) t = 0 and thus 

R _F 3 {r,9) 

~ n 

F 3 (r, 6) X 

P ‘ = —Fr~ 

and furthermore, we can subtract the derivative of (I35al) with respect to r to the derivative of (122al) with respect to p 
and the derivative of (I35al) with respect to 9 to the derivative of (I22bl) with respect to p obtaining (after using (1751) ) 
that (F 3 ) r = {F 3 )e = 0 necessarily, so that F 3 has to be constant and 


(35a) 

(35b) 


Pep 

Pt 


k 

F? 

(36a) 

kX 

(36b) 

772 

r 2 


for some constant k and the arbitrary function F 2 (p + A t) above. 
Integrating (I36al) and (I36bl) . we have 


P (r, 9, p, t) = G {p + At) + F 4 (r, 9) (37) 

where E 4 is an arbitrary function of r and 9 while G is a function of one variable, whose first derivative satisfies 

k 


G' (p + At) = 


(F 2 (p + At )) 2 

while inserting (1251) and (1371) into (I 22 al) and (I22bl) we get the equations 


(F*) r = 


W e 


Aq Be 


{AW + 4A) r 2 rA 2 rB 


(38) 


(39) 


(F) _ r W r rA r rB r 1 
^ i)8 ~ AW + AX 2A + 2B + 2 

and therefore after integration of these last three equations, the dependence of the phase /3 will be explicited. 
To simplify things, we may introduce the function V (r, 9) such that 


F 4 {r,6) :=Vg 


with 


A := 


Be 


2 rVr 


Vw + x 

and such that equation (l39l) is automatically satisfied while equation (l40l) becomes 


v =-Yl. 

v rr — 


1 


Vee 


r 2r 2 r* 

representing an equation for the unknown function V (r, 9). Solutions of equation (|43|) have the form 

V(r,9) :=(^ + Z(r,9) 

in terms of the function £ (r, 9) , and indeed inserting (1U1) into (1751) we obtain the equation 

& , iee 


(40) 


(41) 


(42) 


(43) 


(44) 


A£ = £ 


rr t ^ t ^2 


= o 


(45) 



























iii terms of the Laplacian in the plane r, 6 of the function £ showing that the function £ (r, 8) has to be harmonic in 
the plane of the r, 9 coordinates. Summarizing all the obtained results, we have that the functions 


F 2 (v + Xt) (W + X)*e- 3r ^ 

3 

(46a) 

(rB) 2 

p = G(cp + Xt) + tg 

(46b) 

C = re 2r ^BVW + X 

(46c) 

r p 2r G B 

A:~ ._ 

\JW + X 

(46d) 


+ = (46e) 

satisfy the Dirac equations (I2T1) and the Einstein equations (120fl) , Q20g| ), (1201il) and (I20il) . 

We now solve the remaining Einstein equations: making use of (I44|l . (l45l) and (l46l) . from equations (120al) . (120cll and 
(I20bl) we can derive the second derivatives of B 1 namely 


B rr = A^B ( 4 { ' rB ^ 2 + 8rB - 12 ^ 0d)2 + 4rBB r&0 - WBBg^g 

+5 (rB r ) 2 + 4 B 2 t ee - 2 rBB r + B 2 - 3Bj) 


(47a) 


Bgg = — — (l2 (rB^gf + 8rB 2 £g r g — 4 (B£ggf — 4 rBB r £gg + 4rBBg£ r g 

+3 (rBr) 2 - 4 B 2 (gg + 6rBB r + 3 B 2 - 5 b£) 


(47b) 


B r g = 


ArB 2 £ >r g^gg + 2 B 2 ^ggg + 2 rB r Bg + BBg 
rB 


and thus inserting equations (1371) into (|20c[) . we obtain the expression of the Laplacian 


AW = 


(- kW 3 - 3k\W 2 - 3kX 2 W - kX 3 ) e~ Sr 
r 4 B 2 (W + X) 

[(3 r 4 B 2 + (-12 r%g + 6r 3 ) BB r 


r 4 B 2 (W + X) 

+3 (rBgf + 12r 3 BBg^ r g + (l2r A £ e + 12 (r&gf - 12 r 2 &g + 3r 2 ) B 2 ^ W 2 
+ (6Xr 4 B 2 + ((-24 r 3 &g + 12r 3 ) A - 2r A W r ) BB r + 6A (rBg) 2 

+ (24A r 3 Z r g - 2r 2 Wg) BBg + ((24 r A g e + 24 r 2 $ g - 24 r 2 &g + 6 r 2 ) X 
+4 r 3 W r &g - 4r 3 Wg£ r g - 2r 3 W r ) B 2 ) W + 3X 2 r 4 B 2 
+ ((-12 r 3 (gg + 6 r 3 ) A 2 - 2Ar 4 W r ) BB r + 3A 2 r 2 B 2 e + (l2A 2 r 3 ^e - 2A r 2 Wg) BBg 
+ ((12r 4 & + 12 r 2 $g - 12r 2 ^gg + 3r 2 ) A 2 + (4 r 3 W r ^g - 4 r 3 Wg£ r g - 2 r 3 W r ) X 

+2r 4 W 2 + 2 r 2 Wg-) B 2 ] 

which is a rather heavy expression. Nevertheless, it can be used together with equations (Id7l) into 
and all remaining Einstein equations reduce to the unique final equation 


alt A , 6 rB r £gg 6rBg£ r g 3(rB r ) . 3 rB r 

6 (rCre) + 6 igg - - +- - + —^- G &e 


B 


B 


B 


0 

2 B 2 


(47c) 


(48) 


and (120dl) . 

(49) 
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in terms of B and £ alone. 


2 r£ re 


Be 

B 


rB r 

~B 


-2£ ee + l =0 


and thus splitting into 


rB r 

~B~ 


~ 2 £gg + 1 — 0 


Be 

2 r(i r e + = 0 


the form 


= 0 

(50) 


(51a) 


(51b) 


which can both be easily integrated yielding 


B = 


e- 2r «V 


(52) 


where q is a suitable integration constant. Now, if we set the function W (r, 6) of the form 


W(r,6) 


P(r, 0) 


(53) 


in terms of another arbitrary function P(r,0), and insert (1521) as well as (1551) in the expression of the Laplacian (1481) . 
after some calculations we obtain the Poisson equation 


A P = 



(54) 


for the unknown function P(r, 9). 

And at this point, the complete integration of the coupled system of Dirac and Einstein equations is over. 

For each solution £ of the Laplace equation (l45l) . one gets the function B from equation (l52l) which allows to 
determine the function P(r,6) and thus the function W(r,8) through the Poisson equation (1541) written as 


kr 2 B 2 

r{rP r ) r + Pee = - s — 

for simplicity: inserting the so-found solutions into equations cm, we finally obtain 


(55) 


Vk, 1 
q 3 ^G' (<p + \ij-pi 


(56a) 


P — G + A t) + £g 


(56b) 


A = q 2 VP 

(56c) 

B= q 

r 

(56d) 

ii 

(56e) 

wM =p(ij)- x 

(56f) 
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in terms of the constants k , q, A, the function G and the harmonic function £(r, 9). The corresponding line element is 

dr 2 " 1 


ds 2 = q 4 (2A - A 2 P)dt 2 + 2(1 - A P)d<pdt - Pdq? - e~ 4r ^d9 2 - e~ 4r& 
and the spinor field is given by 


( 57 ) 




y/h e iG(ip+Xt) e i£ e 

^/G'iy+Tt y pF 



( 58 ) 


as the most general exact solution of the spin eigenstates single-handed spinor fields. 

For the sake of completeness, we stress that the determinant of the metric is 

V\detg\ = ( 59 ) 

which will also be used in further applications. 

An interesting feature of solutions (l57l) and (l58l) is that despite both the Riemann tensor and the spinor field be 
non-zero nevertheless all Riemann invariants vanish and so the metric is that of a VSI-spacetime HI and both 
spinor invariants ?/>-!/> and iip-y 5 ^ vanish and so the spinor is a flag-dipole ; in addition to this, with another direct 
check it is possible to see that the momentum density vector defined by V 1 = ’E,~ 2 P~^il>T' l ip is a covariantly constant 
null vector (V 2 = 0 and P{Vj = 0), therefore the metric represents a PP-wave spacetime, according to the definition 
introduced in literature in (7i) (see also [8j). Moreover, the fact that all invariants vanish means they all are regular, 
and this implies that any singularity is removable with a suitable choice of coordinates. 

In the following we are going to study a special case. 


A. The radial case 


A first result that is necessary to stress is a no-go result: the assumption W = 0 amounts to impose non-rotational 
structure to the metric m , and in such circumstance it is immediate to see that from (I56fl) it follows that P is 
constant and then from (15411 we have k = 0 yielding 3 = 0 necessarily, showing that apart from the vanishing spinor, 
no non-trivial solution can exist. This result may appear little in importance, but it shows that we are on the right 
track: the spinor might well be single-handed and a spin eigenstate but still it has a non-trivial rotational structure 
that cannot be compatible with the rotational invariance imposed by the W = 0 constraint; the coupled system of 
Einstein-Dirac field equations is telling in terms of a dynamical solution what was expected in intuitive terms. 

However, this does not mean it is impossible to study the purely radial case. The purely radial case is obtained, 
firstly, by assuming the solution £ of the Laplace equation (1451) of the form in which no elevation angle is present, 
and secondly, by considering also the solution P of the Poisson equation (1541) to have no elevation angle: the first 
assumption gives the possibility to solve the Laplace equation (1451) in general for 


£ = fco + k\ In r 

for some constants ko and k \; this solution together with equation (I56dl) gives 

P r , P e0 ke~ 4kl 


P 


rr \ i o — 2 4 


and this allows us to implement also the second assumption getting 

ke~ 4kl (ln7-) 2 

P = T 0 H-——^-b Ti In r 


2 q 4 


(60) 


(61) 


(62) 


as particular solution, tq and n being suitable integration constants. Inserting equations (1601) and (l62l) into (1561) . we 
have the expression of the spinor 




\/k e iG (v+ At ) 
~<F \JG' {y+W) 



ke 4kl (lnr) 2 

V 


+ T\ In r 


/! 


0 

0 

\0 


( 63 ) 
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and the line element 


ds 2 = - 


q 4 (A 2 ti lnr + A 2 tq — 2A) + 


k (lnr)" A 2 e 


2\2p-4fe! 


df 2 


2q 4 (A Ti In r + Ar 0 — 1) + kX (In r) 2 e 4kl d<pdt 


- ^< 7 4 ti lnr + q 4 r 0 + ^ k (lnr) 2 e 4fcl ^ 


-</ e 


,4 -- 4fcl d6> 2 - g 4 e“ 4fcl 


d(/? 
dr 2 


showing that in the purely radial case, the general solution we found can be fully explicited. 
However, more information is found by choosing a very special tuning for the parameters. 


(64) 


1. The case to = ti = k\ — 0, G(x) = x, k > 0 


A particular instance of the purely radial case m and 
k > 0, where the line element assumes the form 


consists in setting tq = n = k\ = 0, G{x) = x and 


ds 2 = 


2q 4 X- 


k (lnr) A 2 


dr 


2 q 4 — kX (lnr)^ 


dvdt - l%r)_ v - 


(65) 


while the spinor is 


■0 


(^) 4 ^ 

9 2 (In r) ^ 



( 66 ) 


in which we see that In r must be positive, and therefore r > 1 in general. 

This suggest that the singularity of the metric (RjSl) for r = 0 must be apparent, that is depending on the chosen 
coordinates: indeed, by performing the coordinate transformation 


R =: lnr 


(67) 


the line element (16511 assumes the expression 


ds 2 = 


2q 4 X - -kX 2 R 2 


dt 2 + [2<7 4 — kXR 2 ] d(pdt — -kR 2 dip 2 — q 4 d6 2 — q 4 dR 2 


( 68 ) 


which is regular everywhere, as it should have been expected from the fact that no singularity was present for the 
curvature invariants; on the other hand, the spinor field still has a singularity on the sphere of unitary radius r = 1 
or in the new coordinates for R = 0 as clear from 


0 = 


(2 k)i 


,i(ip+\t) _ 


Ri 



(69) 


diverging in the origin, although no singularity is known to be present for the spinorial invariants. It may be possible 
to track the origin of this apparent singularity in the fact that for massless spinors it is impossible to boost in the 
rest frame, and all information about energy and momentum will always be frame-dependent necessarily. 

The computation of the energy of the spinor is done by integrating the energy density over the volume measured 
by the determinant of the metric (1551) and it eventually gives 

E = f ip^il)y/\detg\dV = 4nV2kq 4 f R 
JV J[ 0,oo[ 


( 70 ) 
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where V is the occupied volume: the integral gives no problem at infinity, but in the origin it goes as In and so with 
a positive divergence. This logarithmic divergence recalls the one for ultraviolet regimes in quantum field theory. 

Such divergence in quantum field theory is treated by applying a cut-off and insisting that no calculations be done 
after that, because the theory beyond the cut-off is expected to be different, and in the present context it is possible 
to establish a physical interpretation for this situation: the physics beyond the cut-off is different because there the 
energy density is so large that gravitational effects can no longer be neglected; equivalently, quantum field theory 
ceases to be valid beyond the cut-off because linearity can no longer be used when gravitation is considered. 


IV. CONCLUSION 

In this paper, we have considered the massless neutral spinor field in interaction with its own gravity: as for the 
resulting mutually coupled system of non-linear differential field equations, we have found all of the exact solutions. 

The solutions we have found have been proven to be in the form of PP-waves for the metric tensor and to be a 
flag-dipole for the spinor field: the flag-dipole fields we have found are more specifically the Weyl fermions, which are 
fundamental as the basic chiral components of Dirac fermions; and the property of the metric of being a PP-wave, 
and therefore of describing a VSI-spacetime, implies that such a metric has the radiative corrections that vanish 
identically, and therefore this spacetime after quantization is nothing more than its classical counterpart @. 

Or in a more picturesque language, in this case quantum gravity is nothing more than classical gravity. 

On the other hand the interest for Weyl fermions has very recently received a considerable revival because of the 
discovery of Weyl fermions inside a synthetic metallic crystal called tantalum arsenide, an experimental finding that 
can be used to create massless electrons with no backscattering, therefore improving the efficiency of electronic devices 
such as large-volume single-mode lasers or for general approaches to quantum computing [To| . 

Before applications, it is also from a fundamental point of view that Weyl fermions are important: for instance, it 
is still open the possibility that neutrinos may oscillate because of intrinsic dynamics 0; hence, the persistence of 
the masslessness hypothesis for neutrinos would make them a fundamental example of Weyl fermions. 

If neutrinos were not Weyl fermions because massive in principle, then their masses would most probably assigned 
from Yukawa couplings in terms of the Higgs mechanism after symmetry breaking; but before such a symmetry 
breaking neutrinos as well as all other fermions are still massless: therefore in the earliest stages of cosmic evolution, 
where the universe is supposed to be in a configuration of unbroken electroweak symmetry with all massless fermions 
and dominant gravitational contributions, the particles would be modelled by our solution. 

More in general, our solution is important also from a purely mathematical perspective: because a Dirac fermion 
is a system of two Weyl fermions of opposite helicities with a mixing of the two opposite-helicity states given by the 
mass term then a massive fermion can be thought as a system constituted by two massless fermions plus a sort of 
interaction described by the mass term; consequently, solutions for the massive fermion can be sought in the form of 
exact solutions for the massless fermion plus corrections in small mass regimes. In such a perturbative expansion the 
leading contribution would be the exact solution we have presented here. 

The found solution is therefore important not only because of its technological applications (large-volume single¬ 
mode lasers and quantum computing 0), but also because for the foundations of physics (approaches to quantum 
gravity Q); and of course the solution is important mainly because of the obvious mathematical interest of having a 
non-linear system of differential field equations exactly solved. 

We believe that these results are of high importance for mathematics, as well as for physics in general, consequently 
prompting further research following this direction. 

Extensions are already under investigation. 
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